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Задача оптимального управления компартментних систем 


В даній роботі розглянута задача оптимального керування компартментної системи з анестезіології у 
випадку лінійності та нелінійності з використанням динаміки Міхаєліса - Ментена. 
Ключові слова: оптимальне керування, компартментна система. 


Тре ргобіега ої орійта! сопіго! ої соппрагітепіа! 5узіет ої апезіпезіоїору іп Ше са5е ої Ппеагіїу апі 
попішпеагіїу Бу пзіпє, Фе дупатіся ої Фе МісПасіїз - Мепіеп 15 де8сгібед. 
Кеу ууогаіз: оріта! сопіто!, соппрагітепіа! 5у5кет. 


В данной работе рассмотрена задача оптимального управления компартментного системьі по анесте- 
зиологии в случає линейности и нелинейности с использованиєм динамики Михазлиса - Ментена. 
Ключевьг"е слова: оптимальнос управление, компартментная система. 


Проблеми керування реальними системами у столліття швидкого науково- техніч- 
ного розвитку набувають значну важливість і актуальність. В цьому напрямку викликає 
значний інтерес керування системами з позицій оптимізації певних характеристик си- 
стеми, наприклад, максимізації дальності польоту літального апарату, максимізації при- 
бутку підприємства, мінімізації енергії або витрат, необхідних для досягнення певного 
цільового стану, мінімізації відхилення траєкторії руху системи від заданої та ін. 
Можна навести ще багато подібних задач. Знаходження та дослідження керування, за 
допомогою якого може бути досягнута задана ціль шляхом оптимізації певного крите- 
рію, складає фундаментальну основу теорії керування (1. 

Для постановки задач оптимального керування необхідно, в першу чергу, ви- 
значити цільову функцію оптимізаційного процесу. З цією метою потрібно зробити 
відповідне формулювання задачі у фізичній формі і здійснити переклад цього 
фізичного опису на формальну мову математичних співвідношень |2). 

Задача оптимального керування полягає в тому, щоб знайти таке керування, 
при якому виконуються деякі обмеження 1 мінімізується певний функціонал. Вигляд 
обмежень 1 функціоналів визначаються умовами реальних задач. 
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У даній роботі розглядається задача оптимального керування компартментних 
систем. В таких задачах множина керувань, як правило, наступна: 


С: г) гаки) «б, і, сІХі,, и(1)- вимірна! 


тут а,БЇ, »0 


Нехай припустимо, що стан системи х(/)є К" при заданому керуванні и є(/ 
визначається системою звичайних диференціальних рівнянь: 


ах) 
НО а) 
х(15) 2 Хо 


де /Д:КхК"хК-» К" є неперервною і має неперервні перші частинні похідні 
відносно х та и. Оскільки припускається, що ц(ї) є вимірною та обмеженою, то 


права частина системи (1) є неперервною відносно Х і лише вимірною відносно / 
для фіксованого х. Отже, розв'язки системи (1) є абсолютно неперервними функціями, 
що задовольняють систему (1) майже скрізь. 


Задача оптимального керування містить критерій якості ./ Їм) вигляду: 


Д-р вок у), 


де І - задана дійснозначна функція і ф - неперервно - диференційована 
дійснозначна функція. 

Мета задачі полягає в тому, щоб знайти оптимальне керування и/ є Ї/ таке, щоб 

Днем) 2) 
и 

Після того, як описано модель та визначено критерій якості, в теорії 
оптимального керування ставиться ряд задач: 

-. доведення існування оптимального керування; 

- опис побудови оптимального керування; 

- доведення єдиності оптимального керування; 

-. чисельне обчислення оптимального керування; 

-- дослідження залежності оптимального керування від параметрів моделі. 

Теорема І. Нехай розглядається задача оптимального керування (1), (2) на зада- 
ному інтервалі Ра й Ї: Припустимо, що: 


Г) існує стала М » 0 така, що од «М длявсіхиє Й та І ЗІХІЇ); 

2) І, - напівнеперервна знизу; 

З)множина ФР'- ЩО, у Зу є у з Гху)у 2 Щіхуу є опуклою для 
(1, х) є я Їх ЇМ « м). Тоді існує оптимальне керування и'є 

Теорема 2. Нехай и є ( - оптимальне керування для задачі (1), (2). Тоді існує 
спряжена функція 4: К -» К" така, що х(г,и ),и ,А задовольняють систему: 


ах) Х 
а 90, 6) 
х(5) 2 Ху 
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атакож спряжену систему: 


4) ан ік мій ; 
З --1Ї, Ї, , А Ї, , , 
й Як «би )-А ДО ) 4) 


Ма) я ф'хОТ)), умова ттрансверальності 


де Н - функція Гамільтона - Понтрягіна задається, як: 
Н(іухуи) з Ці,хуи) М Пух). (5) 


Компартментна модель керування 


Розглянемо наступну систему керування: 
«ах 
р - АхОСІ)- /Ох(Т)) є Ви(т) 5 
х(1)є КВ", ицп)є К" 
з відповідною початковою умовою: х(0)-х,. Тут вважаємо, що /:К"-» К" є 
неперервною і /(0)-0. 
Множина керування (/ задається, як: 
Га ша) -Чи, у (0хи ()«10«тхг,, и, (г) -вимірна) 


- кінцевий час керування. 


б 


Метою задачі оптимального керування є експоненціальна стабілізація системи (6). 
Тобто критерій якості в такому випадку має вигляд: 


а Й р'ха)-Кетуз зт 07012 (т) 


И/ -- ваговий коефіцієнт. 

Отже, метою роботи є визначення оптимального керування и є / , ЩО за- 
довольняє умову (2). 

На основі теорема 2 оптимальне керування в задачі (6), (7) існує, оскільки під- 
інтегральний вираз в критерії якості є опуклою функцією, а траєкторія системи на- 
лежить простору 1/7 . Використаємо теорему 2 для отримання необхідних умов опти- 
мальності. 

Для опису побудови оптимального керування використаємо функцію Гамільтона -- 
Понтрягіна вигляду: 


Н-(р'х(І)-Ке")ч опи НА ФАО) я Ох) Ви(т)) 
Мі)є В" 
На основі теореми 2 спряжена система буде мати вигляд: 


СНО ВИН х АР, КРОН СН т 4 (х (т) 
об Хр'х()-Ке""ур-дА (БАчА і (5) 


; - ІН 
На множині у 0 «и, (1)« й маємо п -0,асаме: 
и 


изи 


аН 


-4 заТизА()Ваеб 
ди 


ж 
изи 
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Тоді 


б --ов 


Отже, використавши обмеження зверху та знизу на керування и/ (г), ми прийшли 


до того, що оптимальне керування може бути виражене: 
Т 
: ; В 
и (1) з тіп«, Ся 
| 9) 
Виходячи з теореми І оптимальне керування задачі (6), (7) може бути по- 
будоване в результаті розв'язку такої крайової задачі: 
РД 
. В 
аю Аха)ч РОТ) папа 1,| - дя 


блю В 0) 


Чисельне обчислення оптимального керування 


Методи чисельного розв'язку задач оптимального керування можна класифіку- 
вати як прямі, так і непрямі. Ці методи відрізняються підходами для пошуку розв'язку 
задачі оптимального керування. Непрямі методи намагаються розв'язати крайову задач 
необхідних умов оптимальності. Напротивагу, прямі методи не вимагають безпосеред- 
ньої побудови необхідних умов. Прямі методи не будують спряжену систему, систему 
керування та умов трансверсальності. Вивчаючи оптимальне керування використовують- 
ся обидва підходи. Головним недоліком використання непрямих методів є те, що навіть 
знаючи допустимий стан та керування, немає гарантії, що обчислений розв'язок по- 
кращить відомий. Більше того, непрямий метод потребує початкових наближених 
значень для спряжених змінних, а чисельний розв'язок спряженої системи на практиці є 
слабо обумовленою задачею (31. 

З цієї причини ми використали прямий метод, який дозволяє знайти чисельні 
розв'язки задач, що мають навіть загальніші від (1), (2) постановки. 

Програмна реалізація. Прямий метод чисельного розв'язання задачі оптималь- 
ного керування реалізовано в пакеті /ауа - класів Фуп.Орі. 

Розглядається задача (6) - (11) на основі опису 3-х компартментної моделі з 
анестезіології (4). Опис задачі зроблено за допомогою вхідного текстового файлу, де 
змінні стану системи визначено командою: з/ате хі х2 х3 
змінна керування: сопігої и 
константи: "еа! а1 1 а21 аз1 а!2 а!3 К пи Й рір2р3 
кількість часових вузлів: пойез - 356 
метод розв'язання задачі нелінійного програмування: теїйоа - ауп зар 
метод інтегрування системи диференціальних рівнянь: оає - Пиеп 
точність методу: ерзіїоп - 1.0е-4 

Блок початкових умов: іпійа! сопайіоп: 

хі - 40 

х2 - 0.0 

х3 - 0.0 

У блоці аупатіс едиайоп описується система керування (6) із зазначеними 
параметрами: 

рІ-І 

р2-2.514 
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р3-2.619 

а11-1.6 

а21-0.207 

а31-0.09 

а12-1.901 

а13-1.98 

К-40 

пиз-0. 155 

аа хі --(а!І ча21чаз1) я І чаї 25х2-аї 3 Ях3 чи 

аа х2 - -аї2чх2чаді як І 

ад х3 - -аї3 х3-каЗ І х І 

В блоці іпедиаййу сопзігаїпі описуються обмеження типу нерівності: 

дфе-и 

де иц-0.24 

Критерій якості описується в блоці: со5і /ипсіопа! 

МУ - 0,001 

їпліса! те - 0.0 

па! йте - 56 

ТА (рік рдтх2р31х3)Ч(рІ ЧІ річ 2-р3 їх 3)- 

2Х(рІ ях Ір х2-кр3 3) ЯК Чехри-пи 1) --((К ЯК Зехри-2 пи і)))--0. 5 КИ/ у бу 
240ЕЛ Іедегкі 


192Е-1 


1АДЕ-1 


4792 


000Е0 
000Ео 1Л2Е1 2.24Б1 ЗЗ6Е1 4АВЕЇ 560Е1 


ї 
Рисунок 1 - Поведінка змінної керування дози анестатичного препарату 
пропофолу на заданому часовому проміжку 
400Еї Гедепі 


хі 


о00Ео 1281 2241 З36Е1 4АВЕї 54601 
ї 


Рисунок 2 
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ЗМ 


9А9Е-1 


63ЗЕ- 


ЗЛвЕЧ 


0000 о 
0000 121 


а пет 4АВЕЇ 560Е1 
ї 
Рисунок 3 

На рис. 2 1 рис. 3 відповідно, представлені розв'язки задачі оптимального керування 
при значеннях параметрів та початкових умовах. 

З рис. 1 видно, що на протязі перших 34 хв. слід вводити максимальну кількість 
препарату, а вже на наступному проміжку часу по експоненціальному закону сам пре- 
парат виводиться з організму. 

Оскільки розв'язок лінійної системи можна отримати і в явному вигляді, 
дослідження оптимального керування підсилимо нелінійною задачею з врахуванням 
динаміки Міхаєеліса - Ментена |51. 

Опис задачі, так як і попередньому випадку зроблено за допомогою вхідного 
текстового файлу, де змінні стану системи визначено командою: з/а/е хі х2 х3 
змінна керування: сопігої и 
константи: "еа! а11 а21 аз1 а!2 аї3 КІ К2 КЗ К пи Ирір2р3 

Кількість вузлів, метод та точність розв'язання задачі залишаємо без змін. 
Припустимо, що початкові умови 


хі -40 
х2-00 
х3-00 


Нелінійна система диференціальних рівнянь з врахуванням динаміки Міхаєліса - 
Ментена із зазначеними параметрами, описується в наступному вигляді: 

рі- І 

р2-2.514 

р3-2.619 

а11- 1.6 

а21-0.207 

а31-0.09 

а12-1.901 

а13- 1.98 

КІ-2.7 

к2-1 

к3-0.7 

К-40 

пи-0. 155 

аахі (а! І ча2І чаз1)чх І чаї 29 .МХІАч(КІ -НК2)/к3) на І 3 'ХЗ/х«3-Н(КІ--К2)/К3) ни 

ад х2 - -«аї2чх2-адІ є М/С ЧН(КІ --К2)/К3) 

аа х3- заі3 х3 аз І 'к МОЛ КІ--К2)/К3) 
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Накладаємо обмеження у вигляді нерівностей: 


аф--и 
4д-иц-0.24 


Критерій якості має вигляд: 
Ін(ріїх Ір чх2чр3чх3) рі Ір 2-нр3 їх 3)- 
2ч(рІ ях Ір чх2-кр3 3) ЯК Чехри-пи 0) --((К ЯК Зехр/-2 пи Зі)))--0. 5 КИ ку бу 


240Е1 


1928-41 


ТА 


480Е-2 


ЗЗ1Е9 
000 ЛЛ2Е1 


Шедегкі 


2241 ЗЗ6ЕЇ 4АВЕЇ 5801 
ї 


Рисунок 4 - Оптимальне керування нелінійної системи 
на основі динаміки Міхаєеліса - Ментена 


обо АР 3 


63Е2 


ЗЗЕ2 


249Е-2 


000 й 
000Ео А Різ) 


Іедепа 
хі 


Рисунок 5 


224к1 З.З6ЕЇ 448Е1 540Е1 


Рисунок 6 
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В роботі розглянута задача оптимальної експоненціальної стабілізації. На основі 
принципу максимуму отримано крайову задачу для побудови оптимальної траєкторії 
та керування. 

При побудові чисельного розв'язку застосовується прямий метод, що зводиться 
до задачі нелінійного програмування. 

В якості прикладу розглядається компартментна модель введення анестетич- 
ного препарату. При цьому порівнюється лінійна та нелінійна модель (на основі 
динаміки Міхаєліса - Ментена). З чисельного експерименту видно, що в обох випадках 
приходимо до подібного характеру оптимального керування. При цьому інтегруван- 
ня нелінійної моделі показує на суттєво менші значення концентрації препарату в 
периферичних компартментах, що має суттєве клінічне значення. 
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ВЕ5СМЕ 


У.Р. Мат5епуик, М.М. Сапагунк 
Тие Ргобіет ої Оріїта! Сопіюої! о/ Сотрагітепіаі Зубіет5 


Тре агіїсіє дезбсгібез (ре орійта! сопіто! ргобіет ої сотрагітепіа! 5узіет: 


ах 
а бук ГО) чн Ви) (1) 


х(1)є В", ип)є К" 
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ЗМ 


умі іпісіа! сопаїоп: х(0)-х,. Сопзійег їБаї /:К" -» К"із сопіїпиоц5 апа 


(0-0, 


Сопітої! 5еї 


уже ши) - р оси ()«1,0кх1х1,, ц,(Г) 15 теазигабів| 


і, 15 бпіїе іште сопітої. 


Тре айт ої їБе оріта! сопіго! ргобіет 15 ехропепіїа! зіабійгайоп ої зу5кет (1). 
Собі сгіегіоп їп Різ са5е 15 


РУ РЕ Й р'ха)-Кетуз ит О701 2 


ИЙ із ууеіврійпе Касіог. 

То гергезепі Фе питегіса! 50Їийїоп ої орбіта! сопіго! у/е ц5е Ше фігесі теодй, 
уубісВ 15 шпріетепіеа їп їБе Тауа-сіав5 раскаєе дуп.Орі 

Арріїсабоп ої соптрагітепіа! плоде! ої апезіпейс сретіса! 15 сопзідегед а58 ап 
ехатріє. Іпеаг апа попіїпеаг подеіз аге соппрагед (базед оп Бе дупатіся ої Ше 
Місраєіїя - Мепіеп). Гі Боб са5е5 Ше 5іпаиаг пагиге ої орійта! сопігої м/а5 обіатпей. А 
їБе загае (те іпіергайтя Ше попіїпеаг плоде! 5рому5 зідпійсапіу Їомуег уаїшез Їог Фе 
сопсепігайоп ої сретіса! іп Бе регіррега! сотрагітепія, мбісі 15 ої сппіса! 518 пійсапсе. 


Стаття надійшла до редакції 24.12.2013. 
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